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Abstract 

We present a supersymmetric model of dark energy from Mass Varying Neutrinos which is stable 
against radiative corrections to masses and couplings, and free of dynamical instabilities. This is 
the only such model of dark energy involving fields with significant couplings to any standard model 
particle. We briefly discuss consequences for neutrino oscillations and solar neutrinos. 
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I. INTRODUCTION 



The precision cosmological data acquired over the last decade has changed the nature 
of theoretical cosmology. Models which were previously viable, such as hot or warm dark 
matter, or the open universe model, have been excluded by experiments such as SDSsjlj], 
WMAP[2j, 2dF[3|, and many others fl Q fl. One of the most exciting aspects of the new 
precision cosmology is that we are using gravity to learn a great deal about the dominant 
constituents of the universe. The tremendous results of the High-z Supernovae Search Team 
, as well as the Supernova Cosmology Project (SCP) 8], have given solid evidence that 
the universe is accelerating, potentially due to a new, negative pressure fluid. 

This new fluid, unlike dark matter, does not possess an obvious candidate for any ob- 
servable interactions with standard model particles. In fact, the most popular explanation 
- a cosmological constant - has no dynamics associated with it at all. Slow-roll quintessence 
offers the possibility of detecting an equation of state different from w — — 1, but little else in 
the way of detectable signals. Furthermore, the slow-roll condition requires a field with mass 
on the order of the Hubble scale, 10~ 33 eV. If the quintessence field couples to to standard 
model fields, with couplings of gravitational strength, the typical size of quantum radiative 
corrections to the scalar potential are at least 10 60 times too large, unless exquisitely precise 
cancellations occur. 

Are there theories of dark energy which have naturally sized quantum corrections yet have 
possible non-gravitational signatures (in particular, those arising from couplings to standard 
model fields)? We are willing to assume here the existence of some unknown, perhaps 
nonlocal, physical mechanism which can ensure zero vacuum energy, despite the apparently 
enormous quantum contributions. We still may ask whether, given this assumption, a model 
of dark energy in which the radiative corrections to mass and coupling parameters of the 
model do not have to be finely tuned against the tree level terms. 

In j^l, it was suggested that relic neutrinos could form a negative pressure fluid, and 
yield cosmic acceleration. Such a neutrino fluid could occur if neutrinos interact through a 
new scalar force. The idea is appealing because the neutrino mass scale is comparable to 
that of the dark energy, and the relic neutrinos form a smooth background which we expect 
to permeate the present universe. The result of the scalar force is that the neutrino mass 
depends on neutrino number density, and thus evolves on cosmological timescales. 
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New scalar forces for neutrinos have been considered for quite some time, in the context 



of neutrino oscillations |l0|, , neutrino dark matter [12| , and neutrino clouds 

With dark energy arising from such Mass Varying Neutrinos (MaVaNs), a O(10~ 33 eV) 
mass scalar field is not required. The mass of the scalar field can be as large as O(10~ 4 eV) 
- the scale of the neutrino number density - in order to employ a mean-field approach. 
This mass scale is both comparable to the dark energy scale, and more plausibly stable 



against radiative corrections than the Hubble scale. In ref. 14j| it was argued that the 
local neutrino mass could also depend on the local matter density in this scenario, offering 
the possibility that neutrino oscillation experiments could shed light on the nature of dark 
energy. Much subsequent work has been done on the phenomenology of MaVaNs, including 

Q, 



effects on leptogenesis 

n ■ 



solar neutrinos [18. [19| , the cosmo 
studies of models of various types 



, high energy astrophysical neutrinos [16[ , gamma ra y b ursts |17| , 
ogy of particular forms of the potential |20L |21|, and 
221 ] . Note that other models have attempted to relate 
dark energy and neutrino mass in slow-roll quintessence theories j3, 124I I25I , modifications 
of the energy momentum tensor j3] and a metastable vacuum energy^?!. However, these 
do not rely upon a new, milli-eV scale force between the relic neutrinos, and cannot provide 
the interesting phenomenology of refs , namely, that neutrino mass depends on the 

environment. 

While MaVaN dark energy is an appealing framework, ref. Q, does not present a com- 
plete, radiatively stable model, and questions about the scalar potential and its origins 
remain. In this paper, we provide such a radiatively stable theory of MaVaN dark energy, 
which admits additional interactions with standard model fermions, explains the origin of 
the size of the dark energy in terms of neutrino mass parameters, and can easily mesh with 
a comprehensive effective description of physics below the Planck scale. In section [Hi we 
review the basic features of MaVaN dark energy, including the associated theoretical prob- 
lems. In section we will show that including supersymmetry can not only address the 
quantum corrections, but naturally leads one to consider a new class of dark energy models 
which are analogous to the "hybrid" inflation models 



j . In section IIII A| we describe a 
complete model of neutrino physics, with three additional sterile neutrinos, and discuss the 
consequences for upcoming experiments. In section IIV1 we consider the phenomenological 
implications of this scenario, in particular for solar neutrinos. Finally, in section we 
review the scenario, and discuss future directions. 
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II. DARK ENERGY FROM MASS VARYING NEUTRINOS 

In P] it was proposed that neutrinos interacting via a new scalar force could act as a 
form of dark energy. The idea is most simply understood by promoting the neutrino mass, 
m u , to a dynamical field, which, itself, has an associated scalar potential. The energy of the 
system is then 

V = m u n u + V (m u ). (1) 

Note that m u need not be a canonically normalized field. Assuming the curvature scale of the 
potential is much larger than the Hubble expansion rate, we can use an adiabatic solution 
to the equations of motion where m u tracks the point which instantaneously minimizes the 
total energy, that is, where 

V = n v + V^{m v ) = 0. (2) 

One can then show 

w + 1 = V'(rn„) _ n v 

Vm u fl u + Vt DE ' 

where Qde is the energy density stored in the scalar potential, and Q u is the neutrino energy 
density. 

This can be simply illustrated in a concrete model. We begin by including a scalar field, 
A, which is the canonically normalized field responsible for the dynamical neutrino mass, and 
a fermion N. Because the stored potential energy of A is responsible for the acceleration of 
the universe, we called this field the "acceleron" . N is taken to be a left chiral Weyl field with 
no Standard Model gauge charges, which may be referred to as a "dark", "right-handed", 
or "sterile" neutrino. 

We then consider the interactions, 

CDm D isN + kANN + h.c. + V (A) , (4) 

where k is some Yukawa coupling, and v is a two-component left-chiral Weyl field describing 
some linear combination of the active neutrinos. This Lagrangian is appropriate for energy 
scales well below 100 GeV, as we have integrated out the Higgs field. If kA ^> mo, we can 
remove iV from the low energy effective theory and are left with 

CD^v? + h.c. + V (A). (5) 
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If the relic neutrinos are light mass eigenstates, then we must consider the system at finite 
density, where we have an effective potential for A, 



V cS =n u +V (A). (6) 

Because the first term tends to drive A to larger values, the result of this effective potential 
is to have a neutrino number density dependent value for the A condensate, even if the 
vacuum expectation value is zero. 

A very flat scalar potential Vq, such as a logarithm or a small fractional power, will give 
an equation of state parameter for the acceleron-neutrino fluid which is close to -1, as is 
phenomenologically required for dark energy. 

The framework suffered from two principle theoretical shortcomings. First, a flat potential 
was needed, and both quadratically- and logarithmically divergent radiative corrections were 
too large unless new states appeared in the theory at a scale of order 10~ 2 eV. Second, 
a quadratic potential was not flat enough to give a dark energy with equation of state 
w ~ — 1, and the necessary small fractional power law, or logarithmic, potential had no 
obvious microscopic origin. It has also been recently argued that such theories suffer from 
instabilities at late times j^J. We will return to this point later. 

There are a number of known possibilities to control the size of radiative corrections such 
as supersymmetry, compositeness, extra dimensions, strong near-conformal interactions, and 
pseudo-Goldstone boson type shift symmetries. In this paper, we will focus on supersym- 
metry, which, as a perturbative theory, allows us to consider the theory at arbitrarily high 
temperatures in a theoretically controlled way. 

The question of the origin of the shape of the potential is more difficult. Logarithmic 



31 
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33, 



34j, but only 



potentials do arise in gauge mediated supersymmetry breaking [30, 
at vevs larger than the supersymmetry breaking scale. However, we shall see that a mildly 
different approach will free us of the need for such flat potentials. 



A. Hybrid models 

Many of the theoretical questions associated with dark energy are similar to questions 
that have arisen in the context of inflation. One highly successful framework has been that 
of hybrid inflation In hybrid inflation, the energy density which drives inflation resides 
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in the potential for some field which is in a false minimum due to the large value of another 
light field. A typical hybrid model has a potential 

V = m 2 A 2 + n 2 (Af 2 - v 2 ) 2 + aA 2 M 2 . (7) 

In hybrid inflation, A is referred to as the slow-roll field, while Af is referred to as the 
waterfall field. We use A to make a connection with MaVaN theories. As long as A has 
a large expectation value, Af has a positive mass squared. When the value of A becomes 
smaller than a critical value A c = 2K, 2 v 2 /a, Af becomes tachyonic and will begin to roll. 

In MaVaN theories, rather than being a slowly rolling field, A will have a large expectation 
value due to the presence of relic neutrinos. However, in analogy with the hybrid inflation 
theory, the tachyonic field Af will be stabilized in a false minimum at Af = 0, leaving a 
vacuum energy k 2 v a . So long as the energy density stored in the A condensate is sufficiently 
small, the combined scalar potential will appear as a dark energy with equation of state 
w « —1. 

III. A SUPERSYMMETRIC HYBRID MODEL 

Let us now combine the disparate elements of MaVaNs, supersymmetry, and hybrid 
inflation. Remarkably, in a minimal supersymmetrization of the simplest MaVaN model, 
we immediately arrive at a hybrid model, with the only assumption being the sign of the 
supersymmetry breaking mass squared term for the sterile sneutrino. 

We begin with the simplest Lagrangian for a MaVaN theory, stated in eq. 0] We promote 
v, N and A into chiral superfields £, n and a. The fermionic interactions of eq. 0]are easily 
captured in the following superpotential 

W = nann + m£,£n . (8) 

The resulting scalar potential is 

V = An 2 A 2 Af 2 + n 2 Af 4 + m 2 D Af 2 , (9) 

where Af is the scalar component of n. Here we consider only the light fields in the theory, 
in particular, we set £ to zero as the sneutrino has a large soft mass. Up to a constant, this 
potential bears a remarkable resemblance to eq. [3 other than the sign of the Af mass term. 
However, inclusion of supersymmetry breaking corrections can easily change this sign. 
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Let us consider the supersymmetry breaking quantum corrections to this potential. To 
begin, we should consider radiative corrections above the scale of electroweak symmetry 
breaking. We expect the Dirac mass term to arise from a small Yukawa coupling, W D 
UdHLN . The one loop radiative correction to the M mass is 

5(/,>,)m^-^|^log(A/,), (10) 

where A is the scale at which the soft masses get generated, and m 2 nsy is the average of the 
soft susy breaking higgs and slepton masses squared. Note that for typical values of A and 
m s 2 USJ , the natural size of this correction is of the order m 2 D , and the sign is negative. 

The susy breaking contributions to the M mass squared below the electroweak scale are 
all proportional to K 2 m 2 D , which is in general smaller than mjj, even when log enhanced. 

The acceleron field receives a log enhanced correction to its mass squared of 

Sm 2 A ~ (2m 2 D log(ml/m 2 D ) - 5m 2 N log (A 2 /m 2 D )) (11) 

where fhy is the mass of the active sneutrino, firhj^ is the susy breaking mass of the sterile 
sneutrino, and A is the scale at which Sm^j- is generated. With reasonable values for mo, 
A, and m~, this is roughly ~ 0An 2 (2m 2 D — Sm^). 

Of course, there are also gravity mediated and Planck scale mediated effects whose size 
depend on the gravitino mass and on the degree of sequestering [35] of the a, n sector from 
the supersymmetry breaking sector. We take the calculations of radiative corrections as 
giving a rough lower bound on the natural size for the magnitudes of the soft susy breaking 
masses squared of A and M. 

Altogether, the potential at low energies is 

V = - Cl m 2 D MM* + K 2 N 2 N* 2 + ±k 2 AA*NN* + c 2 K 2 m 2 D AA* + constant, (12) 

where we have assumed the UV physics achieves the appropriate signs, and the c* are 
coefficients of order unity. (Later, we shall see that k is a small parameter, hence we 
ignore terms 0(k 4 ), such as one-loop radiative corrections to the quartic couplings.) 

We can now consider the late time dynamics of this theory, with naturally sized parame- 
ters. We choose the constant in the potential to set the true vacuum energy density to zero. 
The true minimum of this potential is at (A) = and| (J\f) \ = c^rafl/v^K, and the false 
minimum, M = 0, has energy density c 2 mjj/4:K 2 . Since there are no gauge interactions to 
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maintain a large k in the IR, it is natural to assume that k is small. If we wish to identify 
this potential contribution as the present dark energy density, we can thus bound by 
naturalness to be 

m D <i(T 2 - 5 eV. (13) 

Thus, naturalness allows us to identify this as the Dirac mass explaining the solar splitting, 
or some lighter mass, but would involve a one part in 10 3 to 10 4 tuning if we wished to 
identify as the mass associated with the atmospheric splitting. 

If we assume that the relic neutrinos are in the light mass eigenstate, then we can estimate 
the size of the A vev, 



3C(3)T 3 m 
2k 2 k\A\ 

in the non-relativistic case and 



V eS * T + c^mlAA* (14) 



v " - + C2K2m ° AA * (15) 

in the relativistic case. Minimizing the effective potential yields 

kA ~ T/2c 2 /3 (16) 

in the non-relativistic case, and 



kA~ ym D T/12c 2 /2 (17) 

in the relativistic case. We require kA^jTUd, both so the seesaw Majorana mass formula 
works, and, more importantly, so that M has a positive mass-squared. This, in turn, re- 
quires either C2 1 (an unnatural tuning) or mo^zT (meaning the neutrinos are at least 
moderately relativistic). Consequently, naturalness will force the neutrino responsible for 
the dark energy to be the lightest of the three neutrinos, and relativistic. While we might 
enjoy making a direct connection with known neutrino masses, both the scale of dark energy 
and the need for a metastable minimum point us towards the lightest neutrino as the origin 
of dark energy. Note that such theories do not suffer from the potential instabilities of a 
highly non-relativistic neutrino dark energy theory 



A. Multiple neutrinos 

Since there are three active neutrinos, and evidence for at least two mass scales above 
10 _4 eV 2 , we would like to extend this theory to include all three neutrinos and an expla- 
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nation for the atmospheric and solar neutrino oscillations. We restrict ourselves to a single 
acceleron, and three dark fermions iVj. In doing so we should extend the superpotential to 

W = KijariiUj + miiiiii. (18) 

Here we have written the superpotential in the basis where the Dirac mass matrix is diagonal, 
with eigenvalues ordered m 3 > m 2 > m\. The renormalization of the A mass takes place as 
in the previous section, but with additional diagrams possible from the off-diagonal couplings 
between A and the right handed neutrinos. 

The same considerations of naturalness apply here, and we require that the lightest mass 
eigenstate be responsible for maintaining the present false minimum for the A field, and, 
more specifically, that 

4™' > «?X (19) 

where i > j and i > 1 . 

Making this somewhat more quantitative, since T v ~ lCT 4 eV, we must have mi < 10~ 4 eV. 
To achieve the proper vacuum energy, this requires K n ^ 1CT 3 . As a consequence, we expect 
an acceleron mass m A ~ 10~ 7 eV. Taking m 2 > 10~ 2 eV as necessary to explain the solar 
neutrino deficit, and m 3 > 10 1,5 eV as necessary to explain the atmospheric neutrino deficit, 
we have n 2a ~ 10~ 5 and n 3a < 1CT 5 - 5 . 

If the radiative corrections to all the M soft masses are simply proportional to y 2 Dl then all 
masses squared must be negative. If this is the case, then n 3 and n 2 would have rolled off by 
now, achieving vevs rii ~ rui/na (assuming it is stabilized by the quartic piece), resulting in 
an acceleron mass 0(rrii). Since the acceleron Compton wavelength must be longer than the 
cosmic neutrino separation, such a large vev is unacceptable. To reduce the M expectation 
value, one can include terms 

XijknirijUk. (20) 

These need only be large for n 2) 3 in order to prevent a large A mass, and for our purposes, 
we will assume that couplings Xuj to the lightest generation are small. 

Alternatively, it is possible that the corrections to the soft masses squared are not simply 
proportional to y 2 D . One could then have positive masses squared for the heavier two sterile 
sneutrinos, but a negative mass squared for H\. This could arise, for instance, if the heavier 
active sneutrinos are non-degenerate, or in the presence of small corrections form Planck 
scale physics. 
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B. Origin of the Small Couplings 



Our model requires the introduction of several small dimensionless numbers, such as the 
Yukawa couplings which give neutrinos Dirac masses of less than an eV, and the parameters 
Kij. In addition, we have left out several couplings which are not necessarily forbidden by 
symmetries, such as a superpotential term ah u hd, but whose coefficients must be sufficiently 
small. A large number of explanations for such small or absent parameters exist in the 
literature. For instance, one could have a new dimension of size somewhat larger than the 
Planck length, with standard model fields localized to a brane on one end, the acceleron on 
a different brane, and the sterile neutrinos in the bulk, explaining why the acceleron has the 
most suppressed direct couplings to matter and why the couplings of the sterile neutrinos 
have various suppression factors js5 , 



361 ] . In the version of the model with no sterile sneutrino 



vevs, a continuous or discrete lepton number symmetry could suppress unwanted couplings 
such as ah u hd and nhji^. 



C. The Universe from T ~ 1 eV to the Present 



We will assume that the universe is populated by mass eigenstates after BBN. The cir- 
cumstances under which this occurs we will leave for discussion elsewhere 

If the masses- squared for the heavier scalars are positive, then the fluid behaves as a 
traditional hot MaVaN background, but with a quadratic potential, an energy density which 
redshifts far more rapidly than dark energy. 

If the masses-squared for the heavier scalars are negative, then at some critical tempera- 
ture, the A/2,3 fields will become tachyonic and roll from the false minimum. 

The exact vevs and the overall phenomenology is a complicated function of \jk- However, 
we will in general require \\jk <C Xijk^i so that the lightest sneutrino does not acquire a 
large mass when the more massive sneutrinos acquire vevs. For instance, let us consider the 
three-neutrino terms 

W D \(n\ + n\n ?J + n 2 n\ + n\) (21) 

We take the form of the couplings purely for illustration. If we take rh| 3> m|, then we can 
first focus on the early universe behavior of the heaviest sneutrino. 
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M v = 



(23) 



The scalar potential (neglecting couplings to the acceleron) is 

V = -ml | N% | + A 2 1 3A/" 2 2 + + AT 3 2 1 2 + 1 3AT| + 2A2A3 +N%\. (22) 

One can show this is minimized for A2 = =F0.19m3/A, A3 = ±.3m 3 /A. The resulting 
neutrino mass matrix for the heavier two neutrinos is 

/ m 2 \ 

m 3 
777-2 0.6 7773 0.6 7773 

\o 7773 0.67773 1.8 777,3 y 

Under the assumption that 7773 ~ 7773, there are two mostly sterile Majorana neutrinos, with 
masses of 1.8 7773 and 0.6 7773. A neutrino with mass of order O.6777I/7773 is mostly active, but 
with a moderate sterile component. We associate this mass scale with atmospheric neutrino 
oscillations, so the active component is almost entirely ji and r. There is a 1. 7ml/m 3 mass 
neutrino which has very small sterile component, which we associate with solar neutrino 
oscillations. The vacuum energy of the metastable state is A 4 = O.O467773/A 2 . If we take 
7773 ~ 7773 ~ .05eV, the energy density is 0((2 x 10~ 2 ) 4 /A 2 ). 

We now must determine the temperature at which the heavier sneutrino vevs roll to 
their minimum. We will consider the period when all neutrinos are relativistic. Taking the 
acceleron Yukawa matrix to be diagonal, the vev of the acceleron will be determined by its 
mass (which we take by naturalness to be O(Kimi), and the larger of mf/nf for i = 2,3. 
The effective potential is 

V-T^-^ + m^A*. (24) 

If 7773/K 2 > m\j 'K2, then (k s A) 4 ~ K ^I m T 3 ^3 T 2 m\, in which case, the heaviest sneutrino 
would become tachyonic no later than T ~ 7773. On the other hand, if mf//? 2 , > m 4 /^, then 
this can last until T ~ 777 2 , when the intermediate mass neutrinos become non-relativistic. 
In either event, the heavier sterile sneutrino vevs should roll from the metastable minimum 
no later than roughly T ~ 10~ 2 eV (using the solar neutrino data as a guide), but possibly 
much earlier than T ~ .05 eV. 

As a consequence, there will be energy in the oscillating heavier sterile sneutrino field. 
Such scalar field dark matter will cluster like cold dark matter, but with fluctuations 
smoothed out on the neutrino free-streaming length at the time when the field begins to 
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oscillate. It seems difficult to come up with a viable scenario where the entire dark matter 
component of the universe is dominated by such oscillations. It is easy to find parameters 
where the energy density of the oscillating scalar field is much smaller than that of dark 
matter. Even a subdominant dark matter component is very interesting, as it could po- 
tentially be observed or constrained from studies of the matter power spectrum. Such an 
analysis is beyond the scope of this paper. 

Besides the four neutrinos already discussed, the resulting neutrino spectrum will have 
one Majorana MaVaN pair, whose mass continues to evolve on cosmological times. The 
lightest sneutrino remains in its metastable minimum. 

IV. PHENOMENOLOGY 

The neutrino mass phenomenology of this scenario is potentially very rich. With the 
additional singlet states, one must be concerned with active-sterile oscillations, as well as 
active-active, and with the light acceleron and sterile sneutrino fields, novel matter effects 
are possible. We will discuss the two distinct scenarios outlined above: namely, case 1, in 
which there are no sneutrino vevs, has a set of two, heavier, pseudo-Dirac neutrinos with 
a light pair of weakly-Majorana neutrinos, and case 2, which has sneutrino vevs, a heavy 
and light pair of weakly-Majorana neutrinos, and an intermediate mass Majorana neutrino 
whose mass comes from a singlet with a mass of the atmospheric scale. 

In the case 1), we shall see that the splittings of the heavier pseudo-Dirac neutrinos can 
be made sufficiently small that there are no significant consequences, but the splitting of the 
lightest, weakly Majorana neutrinos allows for an MSW enhanced conversion of low energy 
solar neutrinos into sterile states, which must be avoided. We will see that for sufficiently 
small splitting, this need not be a concern, or, in the presence of acceleron-matter couplings, 
the singlet state can be made sufficiently heavy in the sun to prevent level crossings. 

In the case 2), we shall see that the high energy solar neutrinos are not significantly 
changed so long as the singlet neutrinos acquiring masses from sneutrino vevs have masses 
at least 0.03 eV, so as to prevent level crossings with high energy neutrinos. All mixing 
angles involving the two heavier sterile neutrinos are fairly small, but it is possible for the 
new sterile states to affect long-baseline neutrino oscillation experiments. In particular, 
there could be a moderate sterile component in the oscillations of muon neutrinos, with a 
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somewhat shorter oscillation length than that of a — r oscillations. 

The,, have been earlier studies of MaVaNs in the sun. In particular Q studied the 

n 

MaVaN dark energy scenario with the sterile states heavier than an MeV, while [18] studied 
the effects of a particular model in the sun, finding better agreement with data than the 
standard LMA result. However, neither analysis applies to the present discussion. 

There are many other scenarios to consider which we do not discuss further. For instance, 
one could allow direct coupling of the sterile sneutrino fields to ordinary matter. In this case 
the sterile sneutrinos could acquire expectation values which depend on the matter density. 
Since the sterile sneutrino vevs appear in neutrino mass matrices, this could give novel 

n n 

matter effects in neutrino oscillations, as in 14] and discussed in |38J, such as a potential 
explanation of the LSND anomaly 



A. Solar Neutrinos 



In the previous section we saw that a natural theory of dark energy involved three active 
and at least three sterile neutrinos. A natural dark energy mechanism involves one active and 
one sterile neutrino, with Majorana masses in the 10~ 4 eV range. The observed atmospheric 
and solar neutrino oscillations are produced by larger Dirac mass terms. In case 1, none of 
the oscillations into sterile neutrinos affect terrestrial neutrino experiments, as the associated 
oscillation lengths are between 10 5 and 10 8 kilometers for typical neutrino energies. However, 
terrestrial detection of solar neutrinos could be affected by oscillations into sterile neutrinos. 
For solar neutrino phenomena we focus on case 1, as this is potentially more dangerous. 
In this section we discuss solar neutrino phenomena for a particular spectra. We argue 
that the electron neutrinos produced by the B s process observed by the SNO and superK 
experiments are mostly converted to a linear combination of /i and r neutrinos, while lower 
energy solar electron neutrinos are partially converted to sterile neutrinos. 

The terms in an effective Lagrangian involving neutrinos and the acceleron are 

2 

£ D m ai N a Vi + n ab N a N b A + h.c. + ^A 2 (25) 

Here a = 1,2,3 is an index labelling the sterile neutrinos, and % = e, fi,r labels the active 
neutrinos. 

The Dirac type mass matrix m is taken to have eigenvalues « 4 x 10~ 2 eV, rri2 ~ 
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9 x eV, and ni\ ~ 10~ 4 eV. A natural dark energy model may be found where the 
Majorana-type matrix k^ [A) has no eigenvalues larger than of order 10~ 4 eV, has an 
acceleron vev (A) ~ 1CT L0 eV, an acceleron mass m_4 ~ 1CT 7 eV, and in a basis where 
m is diagonal, ku ~ 10~ 3 , k 3ta < and k,2, a ^ 10~ 5 . With three sterile neutrinos, in 

case 1 there will be two pairs of pseudo-Dirac neutrinos with masses of order m 2j 3 and mass 
splittings which are less than of order 10~ 7 eV, and a pair of MaVaN Majorana neutrinos 
with masses and mass splittings of order 10~ 4 eV. The mass rri2 is associated with solar 
neutrino oscillations and m 3 with atmospheric. Because the SNO and SuperK data 40t l41| 
indicate that the more energetic solar neutrinos do not undergo large oscillations into sterile 
states on their way to the earth, we will assume that the mass squared splitting of the second 
pair is less than of order 10~ 12 eV 2 , so that the associated oscillation length is greater than 
1AU for solar neutrinos. We therefore require £22 10~ 8 . 



Due to the MSW effect 
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the high energy neutrinos, with energies above 6 MeV, are 
then mostly adiabatically converted into the active component of the second Dirac pair as 
they propagate from the center of the sun. This component is not exactly a mass eigenstate 
in vacuum, but the splitting between the mass eigenstates is too small to affect neutrino 
oscillations at distances less than many AU. 
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The low energy solar neutrinos are another story. The SAGE [43] and GALLEX 
experiments have detected solar neutrinos with energies from 0.2 MeV, and, when combined 
with other solar neutrino experiments, the results indicate that the survival probability of 
electron neutrinos in the 0.3 MeV energy range should be above 0.47 to be within the 2 
sigma error bars o 



the initial flux 
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the experimental data, assuming the Standard Solar Model is used for 
17] . In our model, however, a large number of these low energy neutrinos 
could potentially oscillate into sterile neutrinos, due to the small mass squared splitting of 
the MaVaN pair and a potential level crossing in the sun from the MSW effect. 

In order to check the acceptability of our scenario for solar neutrinos, we compute the 
effective mass eigenstates and mixing parameters, as a function of the electron density, for 



the following neutrino mass matrix. 



m D 



0.08 5 5 
-0.04 6 40 
0.04 -5 40 



meV, k {A) 



0.1 2 x lO" 5 10~ 4 ^ 
2 x 10~ 5 3 x 10~ 9 10" 9 
^ 10~ 5 10~ 9 10~ 9 



meV (26) 
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The Dirac mass matrix has been chosen to produce a mixing matrix and eigenvalues con- 
sistent with the usual interpretation of neutrino experiments 4^ and the matrix kA was 
chosen to be consistent with our dark energy model, the constraints discussed above, and 
a Frogatt-Nielson j3] explanation of the hierarchy between the couplings of Aq and A^ 2 ,3 
evident in this matrix. 

With these matrices, the heavy pseudo Dirac pair has mass 0.056 eV and mass squared 
splitting of 3.6 x 10~ 13 eV 2 , the second pair has mass 0.0092 eV and mass squared splitting of 
1.7 x 10~ 13 eV 2 and the third pair has masses of 0.00016 eV and 0.00006 eV. In all neutrino 
experiments except solar, the oscillation lengths associated with the sterile neutrinos are too 
long for such oscillations to be observable, and the oscillations among active neutrinos are 
characterized by the mixing matrix 

f 0.84 0.53 -0.09^ 
-0.43 0.56 -0.71 



Vv 



(27) 



y 0.33 -0.63 -0.70 y 



where the rows correspond to flavors e, p, r and columns to masses from lightest to heaviest. 

A full treatment of the effects of the sterile neutrinos on solar neutrino experiments is 
beyond the scope of the paper, but we can make a crude estimate of the effects as follows. 

Flavor evolution of solar neutrinos propagating outward from the core of the sun may be 
computed using the following effective six state Hamiltonian: 

H eR = M 2 /2E + V cc + V nc (28) 

where M is the full six by six Majorana neutrino mass matrix, V cc raises the electron neutrino 
energy by an amount y/2Gpp e , Vnc lowers all the active neutrino e nergies by an amount ^§Pn, 



and p e and p n are the electron and neutron densities respectively |l0j,|42j. In figure 1 we plot 
the eigenvalues of 2EH cS , as a function of x = 2\^2EGpp e , with p n fixed to p e /2. As in the 
usual large angle MSW solution to the solar neutrino problem, at x ~ 10~ 4 eV 2 , a typical 
value for 8 B neutrinos produced in the core of the sun, the third heaviest eigenstate is mostly 
electron neutrino, and this state adiabatically evolves into into a linear combination of e, p, r 
as it exits the sun. The linear combination is not quite a vacuum eigenstate, due to the 
tiny splitting of the pseudo-Dirac pair, but for neutrinos with energies above 0.02 MeV the 
oscillation length is longer than an AU. Thus this state undergoes a nonadiabatic transition 
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to the active linear combination of the third and fourth states, and arrives in roughly equal 
proportions as e, \i and r, as is consistent with the SNO and SuperK experiments. 

The Gallium experiments are sensitive to neutrinos with energies as low as 0.2 MeV, 
which are predominantly produced at x ~ lCT 6 eV 2 . These neutrinos are produced in a 
linear combination of the third and fifth eigenstates of H g r. For large mixing angles and mass 
squared splitting above ~ 10~ 8 eV 2 , the solar evolution is adiabatic [50(. As the neutrinos 
on the third branch exit the sun, they undergo a non adiabatic transition into the active 
superposition of the nearly degenerate third and fourth eigenstates. Since the neutrinos on 
the fifth branch adiabatically evolve into the fifth vacuum eigenstate, the electron neutrino 
survival probability is only 0.23, about 40% the LMA value. 



and Homestake 



These results show the model can easily reproduce the results of the- SNO j^jj], SuperK 



4l| 



5 it ] experiments, but the parameters are likely to rather constrained in order 



for the survival probability of lower energy electron neutrinos to be as high as indicated by 
the Gallium experiments j^, 44, 3|. The basic problem is that in the vacuum, the lightest 



eigenstate is mostly electron neutrino. If this neutrino undergoes an adiabatic level crossing 
in the sun with the lightest sterile neutrino, too few low energy electron neutrinos will 
be seen. One way to avoid this problem is to change the MaVaN parameters so that the 
transition is non adiabatic and suppressed. Either the mixing angle or the mass-squared 
difference can be reduced. In the next section we discuss a more interesting way to achieve 
experimental agreement: addition of an acceleron coupling to matter can raise the mass of 
the sterile neutrino in the sun and avoid the level crossing. 



B. Solar Neutrino Oscillations with matter dependent Acceleron expectation 
value 

In a previous paper, it was noted that a gravitational strength coupling of the acceleron 
to electrons and nucleons could have dramatic consequences for neutrino oscillations jl^ . 
Here we consider the effects on neutrino experiments of the coupling 

C D X N A^-M N , (29) 

Mpi 

where ip^ is the nucleon field, is the nucleon mass, Mp\ is the Planck scale, and the 
above is to be regarded as a term in an effective Lagrangian for energies below the QCD 
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FIG. 1: Plot of the eigenvalues of the matrix 2EH c ff as a function of x = 2\f2Gpp e E u for the 
parameters described in section TlV Al At several values of x, the compositions of the third, fifth 
and sixth eigenvalues are given as (P e , P^ + P T ), where P e ,^,r is the probability that this eigenstate 
will be detected with flavor e, /i, r respectively. The heaviest two eigenstates are a pseudo-Dirac pair 
whose separation cannot be discerned on this plot, and which have negligible electron component. 
The fourth eigenstate is nearly degenerate with the third over most of this region, and is mostly 
sterile. 

confinement scale. The acceleron Compton wavelength is of order 2 meters, and in larger 
objects, the acceleron expectation value may be found by minimizing the effective potential 

V cS = n v E u {A) + V (A) - Xn^P (30) 

Mpi 

where p is the mass density of the object. The second term pushes the acceleron towards 
zero while the first and third terms tend to increase the acceleron expectation value. In 
space, the first term is more important, while for matter at densities above ~ 10~ 6 -^/AjV) 
the second term is more important, and the acceleron expectation value increases. Because 
of the small coupling of the acceleron to the heavier neutrinos, for pAjy < lO 5 ^-, such a 
change in the acceleron expectation value has negligible effect on the two heavier pseudo 
Dirac pairs. The two lighter eigenvalues, however are very sensitive to such a change. 

For 1CT 8 < Xn < 10 3 , the main effect of the coupling (f2HJ) on the evolution of solar 
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FIG. 2: Plot of the eigenvalues of the matrix 2EH c g as a function of y = 2\^2Gpp e , for the 
scenario described in section IIVBI and an energy of 0.3 MeV, which is typical for the Gallium 
experiments. At several values of y, the compositions of the third, fifth and sixth eigenvalues are 
given as (P e ,P^ + P T ), where P e ,^,r is the probability that this eigenstate will be detected with 
flavor e,fj,,r respectively. 

neutrinos is to change the eigenvalues and flavor composition of the two lightest eigenstates 
of the effective hamiltonian 1281 In figures 2-4 we show the eigenvalues of H c g, as a function 
of y = 2y / 2GVp e for three different neutrino energies and p n fixed to p e /2. We have taken 
the acceleron vev to be (A) = (Aq) [1 + 10 r y(MeV/eV 2 )], a crude approximation to the 
expected profile for ~ 1CT 5 . 

We see that the main effect of the inclusion of the term (}2T?|) is raise the mass of the 
lightest sterile state in the sun, and thus to prevent the level crossing at x ~ 10~ 8 eV 2 
between an active and sterile neutrino. At low energies the lightest eigenstate of H e g inside 
the sun is then mainly active. If neutrinos on this branch evolve adiabatically they exit 
the sun as the lightest vacuum eigenstate which is also mostly active. One can thus easily 
achieve a sufficiently high survival probability for low energy electron neutrinos. For optimal 
choice of the coupling A at, one could still have the level crossing with the sterile neutrino 
for intermediate energy neutrinos. This could potentially achieve better agreement with 
experiment than the LMA MSW parameters [52], by substantially suppressing the survival 
of electron neutrinos in the 1 — 6 MeV range, retaining the higher survival probability of 
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FIG. 3: Plot of the eigenvalues of the matrix 2EH e R as a function of y, for the scenario described 
in section HV Bl and an energy of 1 MeV, which is typical for the Homestake experiment. At several 
values of y, the compositions of the third, fifth and sixth eigenvalues are given. 




FIG. 4: Plot of the eigenvalues of the matrix 2EH e ^ as a function of y, for the scenario described 
in section IIVBI and an 8 MeV energy, which is typical for the SuperKamionkande and SNO ex- 
periments. At several values of y, the compositions of the third, fifth and sixth eigenvalues are 
given. 
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lower energy neutrinos, and leaving the survival of higher energy neutrinos unaffected. 

Although we have specialized our discussion to case 1, note that the novel solar neutrino 
phenomenology of the light MaVaN pair of neutrinos is identical for both cases. 



C. New Forces 



For non-zero Xm the acceleron will mediate a new force, which can be searched for in 
tests of gr av ity at distance sea.es shorter than 2 m etets Q Q. In order to eo.pute the 
experimental constraints on X^, the effects of the nonlinear terms in the acceleron equation of 
motion need to be considered, as these can dramatically weaken the constraints |55l . 156 



D. Double Beta Decay 

The consequences of this neutrino mass model for neutrinoless double beta decay are 
simple: there should be no observable neutrinoless decays, even though our neutrinos are 
Majorana. This would be the case even in an inverted hierarchy or degenerate neutrino 
scenario. The explanation is simplest if one works in the flavor basis. All the entries in our 
six by six dimensional neutrino mass matrix are much smaller then the MeV energies typical 
of beta decay. Therefore in computing the double beta decay rate the mass insertion approx- 
imation is excellent. In this approximation the rate is proportional to the m ee component 
of the neutrino mass matrix. This component is zero in the seesaw model. The result that 
neutrinoless double beta decay is unobservable is radically different from the conventional 
seesaw model, where one assumes the sterile neutrinos are much heavier than the beta decay 
energy scale. 

Double beta decay with acceleron or other light scalar emission is allowed but is also too 
suppressed to observe, as only the sterile neutrinos couple directly to light scalars. Thus 
such processes require two insertions of the Dirac mass term, in order to connect the active 
and sterile neutrinos, as well as the Yukawa coupling of the sterile neutrino to the scalars. 
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V. SUMMARY AND CONCLUSIONS 



The MaVaN dark energy model offers the prospect of testable interactions between the 
dark energy sector and neutrinos. However our previous work did not address the issues of 
the origin of the acceleron potential, the degree of fine tuning required, or the stability 
of the dark energy fluid against growth of inhomogeneities. We have shown that in a 
straightforward supersymmetrization of the MaVaN dark energy model, the scalar partner of 
a sterile neutrino acts like the waterfall field of hybrid inflation, whose potential energy can 
drive accelerated expansion of the universe. Supersymmetry can stabilize the required mass 
scales and couplings against fine-tuning, provided that the MaVaN is is lighter than of order 
10~ 4 eV. A 10~ 4 eV MaVaN is sufficiently relativistic to prevent inhomogeneous growth of 
neutrino density. We show that a model with three active and three sterile neutrinos can 
satisfy all experimental constraints while providing a natural theory of dark energy. Possible 
signatures of the model include a reduction of the electron neutrino survival probability of 
low and intermediate energy solar neutrinos, long baseline oscillations of muon and tau 
neutrinos into sterile neutrinos, suppression of neutrinoless double beta decay, and novel 
dependence of neutrino oscillation parameters on the environment. 

As this paper was being completed ref . [58] appeared, which presents a somewhat different 
super symmetric model of MaVaN dark energy. 

Acknowledgements 

The work of A. Nelson was partially supported by the DOE under contract DE-FG03-96- 
ER40956, and a Guggenheim fellowship. N. Weiner was supported by NSF CAREER grant 
PHY-0449818. We thank Alex Friedland, David Kaplan, Cecilia Lunardini and Kathryn 
Zurek for useful conversations. 



[1] K. Abazajian et al. (SDSS), Astron. J. 128, 502 (2004), astro-ph/0403325. 
[2] D. N. Spergel et al. (WMAP), Astrophys. J. Suppl. 148, 175 (2003), astro-ph/0302209. 
[3] W. J. Percival et al. (The 2dFGRS), Mon. Not. Roy. Astron. Soc. 327, 1297 (2001), astro- 
ph/0105252. 



21 



[4] P. de Bernardis et al. (Boomerang), Nature 404, 955 (2000), astro-ph/0004404. 
[5] A. T. Lee et al, Astrophys. J. 561, LI (2001), astro-ph/0104459. 
[6] N. W. Halverson et al., Astrophys. J. 568, 38 (2002), astro-ph/0104489. 
[7] A. G. Riess et al. (Supernova Search Team), Astron. J. 116, 1009 (1998), astro-ph/9805201. 
[8] S. Perlmutter et al. (Supernova Cosmology Project), Astrophys. J. 517, 565 (1999), astro- 
ph/9812133. 

[9] R. Fardon, A. E. Nelson, and N. Weiner, JCAP 0410, 005 (2004), astro-ph/0309800. 
[10] L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 
[11] R. F. Sawyer, Phys. Lett. B448, 174 (1999), hep-ph/9809348. 
[12] M. Kawasaki, H. Murayama, and T. Yanagida, Mod. Phys. Lett. A7, 563 (1992). 
[13] J. Stephenson, G. J., T. Goldman, and B. H. J. McKellar, Int. J. Mod. Phys. A13, 2765 
(1998), hep-ph/9603392. 

[14] D. B. Kaplan, A. E. Nelson, and N. Weiner, Phys. Rev. Lett. 93, 091801 (2004), hep- 
ph/0401099. 

[15] X.-J. Bi, P.-h. Gu, X.-l. Wang, and X.-m. Zhang, Phys. Rev. D69, 113007 (2004), hep- 
ph/0311022. 

[16] P. Q. Hung and H. Pas, Mod. Phys. Lett. A20, 1209 (2005), astro-ph/0311131. 

[17] H. Li, Z.-g. Dai, and X.-m. Zhang, Phys. Rev. D71, 113003 (2005), hep-ph/0411228. 

[18] V. Barger, P. Huber, and D. Marfatia (2005), hep-ph/0502196. 

[19] M. Cirelli, M. C. Gonzalez-Garcia, and C. Pena-Garay, Nucl. Phys. B719, 219 (2005), hep- 
ph/0503028. 

[20] X.-J. Bi, B. Feng, H. Li, and X.-m. Zhang (2004), hep-ph/0412002. 

[21] A. W. Brookfield, C. van de Bruck, D. F. Mota, and D. Tocchini-Valentini (2005), astro- 
ph/0503349. 

[22] R. D. Peccei, Phys. Rev. D71, 023527 (2005), hep-ph/0411137. 
[23] C. T. Hill and G. G. Ross, Nucl. Phys. B311, 253 (1988). 
[24] P. Q. Hung (2000), hep-ph/0010126. 

[25] R. Barbieri, L. J. Hall, S. J. Oliver, and A. Strumia (2005), hep-ph/0505124. 
[26] E. I. Guendelman and A. B. Kaganovich (2004), hep-th/0411188. 
[27] A. Singh, Phys. Rev. D52, 6700 (1995), hep-ph/9412240. 
[28] A. D. Linde, Phys. Rev. D49, 748 (1994), astro-ph/9307002. 

22 



[29] N. Afshordi, M. Zaldarriaga, and K. Kohri (2005), astro-ph/0506663. 

[30] M. Dine and A. E. Nelson, Phys. Rev. D48, 1277 (1993), hep-ph/9303230. 

[31] M. Dine, A. E. Nelson, and Y. Shirman, Phys. Rev. D51, 1362 (1995), hep-ph/9408384. 

[32] M. Dine, A. E. Nelson, Y. Mr, and Y. Shirman, Phys. Rev. D53, 2658 (1996), hep-ph/9507378. 

[33] G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 (1999), hep-ph/9801271. 

[34] A. de Gouvea, T. Moroi, and H. Murayama, Phys. Rev. D56, 1281 (1997), hep-ph/9701244. 

[35] L. Randall and R. Sundrum, Nucl. Phys. B557, 79 (1999), hep-th/9810155. 

[36] N. Arkani-Hamed and M. Schmaltz, Phys. Rev. D61, 033005 (2000), hep-ph/9903417. 

[37] N. Weiner and K. Zurek (2005), hep-ph/0509201. 

[38] K. M. Zurek, JHEP 10, 058 (2004), hep-ph/0405141. 

[39] C. Athanassopoulos et al. (LSND), Phys. Rev. Lett. 81, 1774 (1998), nucl-ex/9709006. 

[40] S. N. Ahmed et al. (SNO) (2003), nucl-ex/0309004. 

[41] M. Miura, Nucl. Phys. Proc. Suppl. Ill, 158 (2002). 

[42] S. P. Mikheev and A. Y. Smirnov, Sov. J. Nucl. Phys. 42, 913 (1985). 

[43] J. N. Abdurashitov et al. (SAGE), Phys. Rev. C60, 055801 (1999), astro-ph/9907113. 

[44] W. Hampel et al. (GALLEX), Phys. Lett. B447, 127 (1999). 

[45] M. Altmann et al. (GNO), Phys. Lett. B490, 16 (2000), hep-ex/0006034. 

[46] J. N. Bahcall (2004), hep-ph/0412068. 

[47] V. Barger, D. Marfatia, and K. Whisnant, Phys. Lett. B617, 78 (2005), hep-ph/0501247. 

[48] M. C. Gonzalez-Garcia and Y. Nir, Rev. Mod. Phys. 75, 345 (2003), hep-ph/0202058. 

[49] C. D. Froggatt and H. B. Nielsen, Nucl. Phys. B147, 277 (1979). 

[50] A. Friedland, Phys. Rev. D64, 013008 (2001), hep-ph/00 10231. 

[51] B. T. Cleveland et al., Astrophys. J. 496, 505 (1998). 

[52] P. C. de Holanda and A. Y. Smirnov (2002), hep-ph/0211264. 

[53] E. G. Adelberger, B. R. Heckel, C. W. Stubbs, and W. F. Rogers, Ann. Rev. Nucl. Part. Sci. 
41, 269 (1991). 

[54] E. G. Adelberger, B. R. Heckel, and A. E. Nelson (2003), hep-ph/0307284. 

[55] J. Khoury and A. Weltman, Phys. Rev. Lett. 93, 171104 (2004), astro-ph/0309300. 

[56] S. S. Gubser and J. Khoury, Phys. Rev. D70, 104001 (2004), hep-ph/0405231. 

[57] B. Feldman and A. Nelson (2005). 

[58] R. Takahashi and M. Tanimoto (2005), hep-ph/0507142. 



23 



